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The differential and total cross sections for kaon pair production in the pp — > ppK + K~ reaction 
have been measured at three beam energies of 2.65, 2.70, and 2.83 GeV using the ANKE magnetic 
spectrometer at the COSY-Jiilich accelerator. These near-threshold data are separated into pairs 
arising from the decay of the 0-meson and the remainder. For the non-0 selection, the ratio of 
the differential cross sections in terms of the K~p and K + p invariant masses is strongly peaked 
towards low masses. This effect can be described quantitatively by using a simple ansatz for the 
K~p final state interaction, where it is seen that the data are sensitive to the magnitude of an 
effective K~p scattering length. When allowance is made for a small number of events where the 
K~ rescatters from the proton, the region is equally well described at all three energies. A very 
similar phenomenon is discovered in the ratio of the cross sections as functions of the K~pp and 
K + pp invariant masses and the identical final state interaction model is also very successful here. 
The world data on the energy dependence of the non-0 total cross section is also reproduced, except 
possibly for the results closest to threshold. 

PACS numbers: 13.60.Le, 14.40.Aq, 25.40.Ep 



I. INTRODUCTION 

There are several mechanisms that can lead to the pro- 
duction of kaon-antikaon pairs in nucleon-nucleon colli- 
sions near threshold. These can be divided mainly into 
two general classes involving (a) the production of a non- 
strange meson that subsequently decays into KK, and 
(b) the associated production of KY* , where the K is 
formed through the decay of the hyperon Y*. The most 
prominent meson that is found almost exclusively in the 
K + K~ /K°K° channels is, of course the 0, and there 
have been measurements of the pp — > pp<f> -^ppK + K~ 
reaction from DISTO Q and COSY-ANKE [|, as well 
as of the pn — > d(f> — > dK + K~ reaction Q. These exper- 
iments show strength also away from the region of the 
peak that could arise from other mesons decaying into 
K+K-. 

The lightest of the non-strange scalar mesons are the 
a${I = 1) and fo(I = 0), which decay mainly into r)Tr and 
7T7T, respectively. Both have masses around 980MeV/c 2 
and widths of the order of 50-100 MeV/c 2 0J. Since 
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the central mass values fall almost exactly at the KK 
threshold, the strong coupling to this channel distorts 
significantly the upper parts of the mass spectra [f| • It 
has even been suggested that both resonances might be 
mainly molecular in nature [||. Since the K + K~ system 
is a mixture of isospin 1 = and I = 1 , one can consider 
whether the non-0 events observed in pp — > ppK + K~ 
might be linked to the production and decay of the ao 
and/or the /o mesons. The 1=1 channel is isolated 
cleanly by looking at the pp — ► dK + K° reaction Q. 

On the other hand, we can also expect kaon pairs to 
be created in the pp — > K + pY* reaction, where the hy- 
peron decays through Y* — > K~p. There are, of course, 
several excited hyperons that could contribute to such 
a process. Of particular interest for low energy produc- 
tion are the £(1385) and the A(1405). Though nom- 
inally lying below the sum of the K~ and p masses, 
their large widths (~ 50MeV/c 2 ) ensure that they over- 
lap the K~p threshold Results have recently been 
presented on the production of these states in the pp — > 
if+p£ (1385) and pp — > K + pk{lA0h) reactions, where 
the hyperons were detected in the A7r° and £°7r chan- 
nels, respectively Q. 

In addition to being involved in a direct production 
mechanism, the A(1405) and the £(1385) might also lead 
to a significant K~p final state interaction (fsi) in the 
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PP ~ * ppK + K~ reaction. Evidence that the K~p fsi is 
indeed important is to be found in the measurement of 
the pp — > ppK + K~ cross section at energies below the 
threshold Q. It was shown that the invariant masses of 
the K~p final system were on average lower than those 
of the K + p pairs and it was suggested that this is due to 
the strong attraction in the K~p system. 

It is therefore an interesting but open question whether 
the pp — > ppK + K~ data from the non-0 region should 
be viewed, in first approximation, as scalar resonance 
production via pp ppao/fo or hypcron production 
through pp -» !f+pA(i405)/E (1385). It is one aim of 
the present paper to study this competition by using the 
data obtained when investigating 0-meson production in 
proton-proton collisions at ANKE Q. 

A program to investigate the production of kaon pairs 
in nucleon-nucleon collisions has been initiated at the 
COSY storage ring of the Forschungszentrum Jiilich. Al- 
though the methods of the measurement and analysis 
of the pv—> ppK + K~ reaction have already been given 
in Ref. [2j , for the sake of clarity the main points are re- 
peated in Sect HIl Due to the finite width of the 0, distin- 
guishing between events corresponding to the production 
of this meson and the rest requires a detailed modelling 
of the K + K~ spectrum, as described in Sect. Mil When 
more than one pair of particles interact strongly in the 
final state, there is no reliable prescription to evaluate a 
corresponding enhancement factor. Our approach to this 
problem in terms of an effective K~p scattering length a 
is described in Sect. HVi with a heuristic justification be- 
ing given in Appendix [A"l 

The relative distributions of the K~p and K + p in- 
variant masses are shown and discussed in Sect. |V| The 
ratio of the differential cross sections as functions of 
the K~p and K + p invariant masses indicates that the 
K~p attraction is very strong. This confirms the ear- 
lier COSY-11 findings [§] but with higher statistics over 
much wider ranges of energy and mass. The events that 
come from the non-0 region in the 2.65 GeV data al- 
low us to fix the parameter a necessary to describe the 
K~p/ K + p ratio within our final state interaction model. 
The data are mainly sensitive to the magnitude \a\ and 
the value obtained is not dissimilar to that required to 
describe free K~p scattering. 

The approach reproduces also our results at the other 
two energies as well as the COSY-11 data. If allowance 
is made for the small possible contribution to this ra- 
tio from events, of the order of 20%, the two-particle 
ratio is also well described for all K + K~ mass inter- 
vals. In addition, we find a completely analogous effect 
in the three-body Kpp system, with the K~pp/ K + pp 
ratio being strongly biased towards the lowest values of 
the invariant mass. We show in Sect. IVII that these dis- 
tributions are also predicted quantitatively by the same 
K~p final state interaction model with the same value of 
the parameter a. 

The total cross sections for and non-0 production 
are presented in Sect. IVIII The energy dependence of 



the non-0 total cross section is strongly influenced by 
both the pp and K~p fsi but, even after including these 
effects, the lowest energy data are not well reproduced. 
This leads to the suggestion that the effect might be con- 
nected with the two or three points in the K + K~ spec- 
trum at very low invariant masses, which lie consistently 
above the simulations. Our conclusions and views on the 
outlook are presented in Sect. I Villi 



II. EXPERIMENT AND DATA ANALYSIS 

The study of kaon pair production has been performed 
using data taken with the ANKE magnetic spectrome- 
ter [ijj]- This is placed at an internal target station of 
COSY, which is the COoler SYnchrotron of the Research 
Center Jiilich [ll|. We have already reported an initial 
analysis of the pp — ► ppK + K~ data, from which values 
of the pp — > pp(j> cross sections were extracted 0] . 

The measurements were performed at three proton ki- 
netic energies of T p = 2.65 GeV, 2.70 GeV and 2.83 GeV, 
corresponding to excess energies of s = 51 MeV, 67MeV 
and 108 MeV with respect to the ppK + K~ threshold. 
A dense hydrogen cluster jet target [l^] provided areal 
densities of ~ 5 x 10 14 cm" 2 and the average luminosity 
during the experiment was typically 2.3 x 10 31 cm _2 s -1 . 
The detection system of the three-dipole magnetic spec- 
trometer ANKE registers simultaneously slow positively 
and negatively charged particles, with fast positive par- 
ticles being measured in the forward system [Hi EH, EH • 

Kaons produced in the pp—>ppK + K~ reaction at 2.65 
and 2.70 GeV have laboratory momenta between 0.2 and 
l.OGeV/c and this full range was covered by the detec- 
tion system for both K + and K~ . For the 2.83 GeV run, 
the K + were detected only up to 0.6GeV/c, which re- 
duced the geometrical acceptance. The momentum range 
of the protons emerging from the reaction is between 0.8 
and 1.6GeV/c, for which the forward detector has full 
acceptance. 

The momenta of both the kaons and protons were de- 
termined from the track information furnished by sets of 
multiwire proportional chambers. The tracking efficiency 
for kaons fell smoothly from 98% to 93% as the momen- 
tum increased, whereas that for protons was about 80% 
over the whole momentum range. The efficiency map, in- 
cluding its dependence on the particle momenta and posi- 
tion, was generated for the track reconstruction and used 
in the analysis. The data-taking efficiency was around 
90% and, taken in combination with the tracking effi- 
ciency, resulted in a total detector efficiency of about 
60%. 

In order not to reduce the geometric acceptance un- 
necessarily, only one forward-going proton was mea- 
sured in coincidence with the charged kaon pair. The 
pp^>ppK + K~ reaction was then identified by requiring 
that the missing mass in the reaction was consistent with 
that of the non-observed proton. As a first step, positive 
kaons were selected through a procedure, described in 
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detail in Ref. [H, that used the time of flight (TOF) be- 
tween START and STOP scintillation counters of a dedi- 
cated K + detection system. In the second stage, both the 
K~ and forward-going proton were identified from the 
time-of-flight differences between the STOP counters in 
the negative and forward detector systems with respect 
to the STOP counter in the positive system that was hit 
by the K + . By comparing this with the time difference 
expected for a kaon pair based upon their measured mo- 
menta, a clear separation is achieved, with good events 
lying in an island where the two time differences are con- 
sistent [TJ. These two TOF selections, as well as that 
for the K + , were carried out within ± 3er bands. 
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FIG. 1: Missing-mass distribution of pp^K + K~pX events 
at T p =2.65 GeV. The hatching indicates the range used for 
the selection of protons. The side-band events scaled by the 
solid line are indicated by the dotted shading. 



The spectrum of missing masses with respect to the 
K + K~~p system, shown in Fig.Q]for the 2.65 GeV data, is 
dominated by the proton peak. The corrected numbers of 
K + K~pp events obtained in this way are around 3100 (at 
T p = 2.65 GeV), 1300 (2.70 GeV), and 650 (2.83 GeV). 
The estimated background inside the proton cut window 
is 5%, 12%, and 18% at these three energies, respec- 
tively. This has been subtracted using the side-band 
events scaled by the solid line shown in Fig. [TJ The un- 
certainty resulting from this procedure is estimated to 
be less than 3% for the total cross sections and negligi- 
ble compared to the statistical errors for the differential 
distributions. 

After the identification of the pp — > ppK + K~ events, 
the data were binned in intervals of the K + K~ invariant 
mass, with the 2.65 GeV results being presented in Fig. [5] 
This spectrum shows a very prominent <p peak positioned 
close to the nominal mass of the meson. In addition there 
is a non-</> contribution that decreases steadily from low 
to high invariant masses. This behavior, which is mainly 
a reflection of the ANKE acceptance for non-resonant 
ppK + K~ production, has to be simulated in order to 
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2: The K + K~ invariant -mass distribution for the 
pp — > ppK + K~ reaction at 2.65 GeV. The prominent <f> peak 
sits on a steadily falling background corresponding to the 
ANKE acceptance for the four-particle final state. The sim- 
ulation of these two contributions and their sum are shown, 
respectively, by the dotted, dashed, and solid histograms. 



The acceptance has been determined through a Monte- 
Carlo simulation based upon GEANT4 where the 
geometrical acceptance, resolution, detector efficiency, 
and kaon decay probability were taken into account. In 
the initial step, the (^-production distributions previously 
published [|J were taken as the basis of the simulation, 
whereas four-body phase space was used for the non-(/> 
component of Fig. O There are seven degrees of freedom 
for a four-body final state produced in an unpolarized 
reaction. Thus for each beam energy, seven independent 
cm. distributions were generated (five angular distribu- 
tions, the relative momentum of final protons, and the 
K + K~ invariant mass). These were then divided into 
two groups depending on the value of the K + K~ invari- 
ant mass, i.e., </>-rich and <^>-poor regions. All distribu- 
tions were jointly fitted to the experimental data and the 
relative contribution of <p an d non-<^> production evalu- 
ated for the determination of the acceptances. 

In the 0-poor region, most of the cm. distributions 
are fairly uniform, though the polar angle of the K + K~ 
system relative to the beam axis in the overall cm. sys- 
tem shows deviations. Fitting this with the function 
1 + a cos 2 6* gives anisotropics of a = 0.5 ± 0.2 and 
a = 1.5 ± 0.5 at 2.65 GeV and the two higher energies, 
respectively. The inclusion of this factor increases the 
total acceptance for non-<^ production by about 10%. In 
addition, as will be seen in Sects. |V] and IVI| the A' ± p 
and K^pp invariant masses deviate strongly from phase 
space. These deviations have been taken into account 
in the simulations, using iteratively the theoretical ap- 
proach discussed in Sect. IIVI in order to converge on the 
acceptance-corrected distributions. The simulations of 
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the an d non-0 components of the uncorrected K + K~ 
spectrum at 2.65 GeV are shown in Fig. [5] together with 
their sum. 

The estimated total acceptance for kaon-pair pro- 
duction is 1.1% (2.65 GeV), 0.7% (2.70 GeV) and 0.1% 
(2.83 GeV), and the associated systematic uncertainties 
in the total cross section are 5%, 7% and 12%, respec- 
tively. The shapes of the differential cross sections, es- 
pecially those for the Kp and Kpp invariant masses, are 
very stable to this iterative procedure. As a consequence, 
the uncertainty that this introduces in the differential 
cross sections is negligible compared to the statistical er- 
rors. 

The luminosity required to evaluate absolute cross sec- 
tions was determined by measuring in parallel proton- 
proton elastic scattering for laboratory angles between 
5.0° and 8.5°. For this purpose the momentum of a 
forward-going proton was determined using the ANKE 
forward detector. After a missing-mass analysis, the 
(small) background underneath the proton peak was sub- 
tracted and the remaining events were corrected for effi- 
ciency and acceptance. 



III. THE K+K- INVARIANT MASS 
DISTRIBUTION 

Figure [3] shows the pp — > ppK + K~ differential cross 
sections in terms of the K + K~ invariant masses as ob- 
tained at our three energies. These have been corrected 
for acceptance, deadtime, etc., and normalized on the 
basis of the proton-proton elastic scattering data, as dis- 
cussed in Sect. IVIII The resulting non-0 contribution 
is much flatter than that illustrated in Fig. O In order 
to ensure a clean separation of the from the rest, it is 
necessary to model the shapes as reliably as possible. 

In our earlier work Q , the non-0 contribution was as- 
sumed to be described by a four-body phase space. The 
was taken to have a Breit-Wigner form with a constant 
width of r — 4.26MeV/c 2 4], which was then convo- 
luted with an experimental resolution of a = 1 MeV/c 2 . 
Only for the component was the spectrum modified 
to include effects arising from the final state interaction 
between the emerging protons. These three assumptions 
have now to be reconsidered. 

Since the is a p-wave resonance that decays strongly 
into KK, its width depends upon the K + K~ and K°K° 
relative momenta. Due to the narrowness of the state, 
the peak itself is little affected by this refinement but 
it does suppress slightly the low mass tail of the reso- 
nance. Although for low K + K~ effective masses there is 
no clear evidence for any pp fsi, only minor effects are to 
be expected there because of the large amount of energy 
that is then available to excite the pp system. Neverthe- 
less, the fsi cannot be neglected at higher K + K~ masses, 
where its inclusion distorts slightly the four-body phase 
space spectrum. Finally, as will be seen in Sect. [Vj it is 
clear that there is a large final state interaction between 
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FIG. 3: Differential cross section for the pp ^ ppK + K 
reaction with respect to the K + K~ invariant mass at (a) 
2.65 GeV, (b) 2.70 GeV, and (c) 2.83 GeV. The dotted his- 
togram and the dashed curve represent the simulations for 
the <j> and non-0 contributions to the spectrum, with their 
sum being shown by the solid histogram. For all three ener- 
gies the simulations fail to reproduce the lowest mass points. 



the K~ and one or both protons and this also modifies 
the shape of the non-0 component, enriching the central 
part of the K + K~ spectrum. By chance, all three effects 
go in the same direction. In the fits to Fig. [3J the frac- 
tions to be associated with production are diminished 
somewhat compared to those reported in our previous 
work 0]. 

It must be noted that the lowest mass points lie above 
the fitted curves for all three data sets of Fig. [3l These 
deviations from the simulations are significant because, 
as seen from the comparison of Figs.[2]and[3^,, the ANKE 
acceptance is largest in this region. Although ignored in 
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the fitting process, they were taken into account in the 
evaluation of the non-c/i contribution to the total cross 
section, contributing 4% (2.65 GeV), 6% (2.70 GeV), and 
10% (2.83 GeV). 



IV. FINAL STATE INTERACTION 
ASSUMPTIONS 

We wish to describe the pp — ► ppK + K~ data at all en- 
ergies in terms of final state interactions between the two 
protons and between the K~ and either one or both of 
the protons. There is no reliable way of carrying out such 
a program without solving multibody equations, which is 
well outside the scope of the present paper. As explained 
in Appendix [Al we make the ad hoc assumption that the 
overall enhancement factor is the product of the enhance- 
ments in the pp and two K~p systems, all evaluated at 
the appropriate relative momenta q: 

F = F pp (q pp ) x Fx p (qK Pl ) x FK p (qK P2 ) ■ (1) 

Similar approaches are to be found elsewhere in the liter- 
ature as, for example, in the description of the pprj final 
state at low energies jl7j ] . 

The K~p enhancement factor is taken in the scattering 
length approximation 

1 



F Kp {q) 



(2) 



1 — iqa 

where some caution should be exercised in the interpreta- 
tion of the complex parameter a as an effective scattering 
length. In principle a could have an energy dependence, 
which would correspond to the introduction of an effec- 
tive range term. However, it is hard from our data to 
determine the values of any extra free parameters. 
A proton-proton enhancement factor of the form 

\F PP (l)\ 2 = 4^4 ( 3 ) 
q z + a 1 

is assumed [l8[ , where the position of the x So virtual state 
is well fixed at a = 0.1 fru -1 . The final state interaction 
should not influence higher partial waves and, in order 
to reduce its effects above about lOMeV, we take (5 — 
Q.5fm . The relative momentum of the two protons in 
the pp rest frame, which is very sensitive to the form of 
the pp fsi, is then reproduced reasonably, as seen from 
Fig. 1 

It is important to note that the form of neither the pp 
nor the K~p enhancement factor is valid for an excitation 
energy below about 1 MeV due to the explicit Coulomb 
force. However, these regions represent only very small 
fractions of the allowed phase space. 



V. THE KAON-PROTON INVARIANT-MASS 
DISTRIBUTION 

In data taken with the COSY-11 spectrometer at 
e = 28 MeV, i.e. at an energy where there is very lit- 
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FIG. 4: Relative momentum of the final protons from the 
pp — > ppK + K~ reaction in the pp rest frame. The experimen- 
tal points at 2.65 GeV are compared to the simulation, using 
the standard parameters, which is shown as a histogram. 



tie <f> production, it was found that the K~p and K + p 
invariant-mass M(Kp) distributions were markedly dif- 
ferent @ . The ratio of the acceptance-corrected distri- 
butions 



R 



Kp 



dafdM(K-p) 
da/dM(K+p) 



showed a very strong preference for low values of M(Kp). 

Figure [5] shows the Kp spectra obtained at our low- 
est energy. In order to consider a situation similar to 
that of COSY-11, the data have been selected such 
that only events with K + K~ invariant masses less than 
1.01 GeV/c 2 are retained. The simulation that takes into 
account only the pp fsi fails to describe the data. Further- 
more, the ratio of the K ~~p to K + p mass distributions is 
far from constant, as it would be if there were only the pp 
fsi. The peaking of Rkp to the lowest invariant masses 
confirms the COSY-11 observation Q but with higher 
statistics and over a wider range of Kp masses. 

When the K~p final-state-interaction factors are also 
introduced into the four-body phase space simulation of 
the non-^> contribution to the pp — > ppK + K~ reaction, 
the individual da/dM{K+p) and da/dM(K~p) distri- 
butions of Fig. [5l as well as their ratio Rkp, are all well 
described. The best fit to these data is achieved with 
an effective K~p scattering length of a — (0 + 1.5i)fm 
and this value of a will be retained for all the subsequent 
simulations. Nevertheless, the uncertainties in the real 
and imaginary parts are quite large and strongly corre- 
lated such that | a | is determined much better than its 
phase, as illustrated by the \ 2 contour plot of Fig . [5] 
The free K~p I — and 1=1 values from Ref. [19| are 
also shown there along with their one-cr limits. From 
this it is seen that there is an overlap with either isospin 
scattering length. One has also to bear in mind that our 
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FIG. 5: Differential cross sections for the pp — > ppK + K~ re- 
action with respect to the invariant mass of (a) the K + p and 
(b) the K~p systems. In panel (c) their ratio is shown on a 
logarithmic scale. These 2.65 GeV data have been selected to 
come from the region M{K+K~) < 1.01 GeV/c 2 , where there 
is only a very small (j> contribution. The dashed histograms 
represent the results of four-body phase-space simulations, 
where the only distortion is that coming from the pp fsi. The 
solid histograms include also the K~p fsi. 



data represent averages over some energy range and that 
a constant scattering length is necessarily an oversimpli- 
fication. 

The situation for higher K + K~ masses is complicated 
by the contribution from the <fi meson. At our beam en- 
ergies, the average separation of the K~ (produced in 
the tfi decay) from one of the final protons is about 7 fm. 
This large distance reduces considerably the effects of 
the K~p fsi. We would therefore expect that deviations 
of Rkp from unity should arise primarily from the non-0 
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FIG. 6: Contour plot corresponding to the one-cr level un- 
certainty of the fits of the complex parameter a of the fsi 
ansatz to the data of Fig. [5] Our best fit is obtained with 
a — 1.5i fm (marked with a cross). Also shown are the best fit 
values (points) and corresponding contours (dashed lines) for 
the isoscalar (ao) and isovector (ai) K~p scattering lengths 
derived from the study of free 7f~-nucleon scattering [ljj]. 



contributions. The modelling of the data thus depends 
in the first instance upon the determination of the frac- 
tion of 4> production in the overall pp — ► ppK + K~ spec- 
tra presented in Fig. [3] However, classical simulations 
suggest that more than 15% of the <f> decays lead to at 
least one K~p pair with a separation of less than 1 fm, 
for which the fsi is not negligible. 

To investigate this point further, the 2.65 GeV data 
have been divided into four regions of K + K~ invari- 
ant mass, viz. M(K + K~) < 0.995 GeV/c 2 , 0.995 < 
M{K+R-) < 1.003 GeV/c 2 , 1.003 < M{K + K ) < 
1.010 GeV/c 2 , and 1.010GeV/c 2 < M(K + K~). The first 
two intervals contain almost no <fi events, whereas the last 
corresponds mainly to <f> production. The correspond- 
ing Rk p distributions are shown in Fig. [7l where they 
are compared to the same simulation as that shown in 
Fig. [5l The non-0 region is very well described but, for 
high K + K~ invariant masses, the asymmetry seems to 
be diluted too much by kaons from the cf> decay. The 
agreement is improved by assuming that 20% of the <fi 
events are influenced by the K~p fsi. The same effect 
could be achieved by artificially reducing the fraction of 
(j> mesons in the peak of Fig. [3^ by 10%, but this would 
be hard to justify. 

If the strong K~p final state interaction is responsible 
for the distortion in the K~p/ K + p ratio seen in Fig. [7j 
similar effects should be seen for the other excess ener- 
gies. For our e — 66.6 and 108 MeV results shown in 
Fig. the statistics are lower and we merely divide each 
data set into two, comprising 0-rich and 0-poor regions. 
Without any change in the value of the K~p effective 
scattering length, all the distributions are well described, 
especially if it is assumed that 20% of the <f> events lead 
to K~p fsi effects. 

Data on the K~p/ K + p ratio were not published by 
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FIG. 8: Values of the Rk p ratio for (a) the 0-poor region 
and (b) the 0-rich at 2.70 GeV with (c) and (d) showing re- 
spectively the analogous results at 2.83 GeV. The predictions 
obtained using the standard parameters are shown with no fsi 
effect for the K~ from the decay of the cf> (solid curve) and 
20% (dashed). 



shown by the 28 MeV data is convincing. Unfortunately, 
in the publication [9( the values of the ratio at 10 MeV 
were scaled upwards to make them lie closer to those 
taken at 28 MeV. Since the numbers of K~p and K + p 
events are the same, the weighted average of the ratio 
must be unity and this allows one to make an informed 
guess regarding the value of the scaling that has to be 
reintroduced. No account has been taken of the COSY- 
11 invariant-mass resolution, which will be of greater sig- 
nificance for the 10 MeV data. 



FIG. 7: Values of the Rk p ratio for the pp —> ppK K~ 
reaction at 2.65 GeV extracted for four different selec- 
tions of the K + K~ invariant mass, viz. (a) M(K + K~) < 
0.995 GeV/c 2 , (b) 0.995 < M(K+K~) < 1.003 GeV/c 2 , (c) 
1.003 < M(K + K~) < 1.010 GeV/c 2 , and (d) 1.010GeV/c 2 < 
M(K + K~). The histograms are predictions of the model de- 
scribed in the text, where the K~ mesons from non-<^> events 
are subject to a K~p final state interaction. In the high mass 
interval, the effects of assuming a 20% probability of the K~ 
from the (j> decay being influenced by a fsi is considered and 
this leads to the dashed histograms. 



the DISTO collaboration [l[. However, the same good 
description of Rk p is found for the COSY-11 results of 
Fig. [9] where, in view of the low excess energies, there is 
no ambiguity regarding the <f> contribution. Given that 
the systematic uncertainties in this experiment are inde- 
pendent of those at ANKE, the quantitative agreement 
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FIG. 9: The R Kp ratio for the COSY-11 pp -> ppK + K~ data 
at (a) e = 10 MeV (upper panel) and (b) e = 28 MeV. The 
small contribution from the low mass tail of the <f> has been 
neglected in the 28 MeV simulation. Since the 10 MeV results 
were arbitrarily scaled in the publication they have here 
been divided by a factor of three before being compared to 
the predictions of the fsi model, using the same values of the 
parameters as for the ANKE data. 



VI. THE KAON-PROTON-PROTON 
INVARIANT-MASS DISTRIBUTION 

We have argued in Sect. IIVI and Appendix [5] that it 
is natural to take the final state enhancement to be the 
product of the individual factors corresponding to the pp 
and two K~p interactions, which allows the kaon to be 
attracted simultaneously to both protons. This in turn 
suggests that there might also be an enhancement in the 
three-body system at low K~pp invariant masses. This 
is indeed the case. The size of the effect at 2.65 GeV is 
illustrated in Fig. [10] for the regions with M(K + K~ ) ^ 
1.01GeV/c 2 . 

There is a striking preference for low Kpp masses 
in both K + K~ intervals and this behavior is repro- 
duced quantitatively by our simulation, which uses the 
same K~p effective scattering length a that was chosen 
to describe the two-body Rkp ratio. It provides some 
empirical justification for our choice of the final-state- 
interaction enhancement factor as the product of the in- 
dividual two-body enhancements. This favors configura- 
tions where the K~ is attracted simultaneously to both 
protons. It must, however, be stressed that the strong 



mass dependence which is apparent in Fig. [10] does not 
necessarily imply that the K~ forms a bound state with 
the two protons [H [M [H . 

Though the error bars are considerably larger at the 
two higher energies, as seen in Fig. 1111 the simulation 
is consistent with the results. In all these three-particle 
mass distributions the introduction of a small amount of 
fsi for the products of the c/> decay seems to be of less 
importance than for the Kp spectra. 

VII. TOTAL CROSS SECTIONS 

The total cross section results for both <f> and non- 
<f> production are given in Table U for the three energies. 
The systematic errors given here arise from the quadratic 
sum of the uncertainties in the background subtraction, 
acceptance correction, tracking efficiency correction for 
kaons (3%) and the analysis of the pp-elastic data used 
for the normalization. Unfortunately, the pp database is 
very limited in our region of angle and energy. This is 
recognized in the recent SP07 update to the SAID pro- 
gram where it is stated that our solution should be consid- 
ered at best qualitative between 2. 5 and 3 GeV (23j . The 
raw SP07 solution generally overestimates the values of 
those pp differential cross sections that are available in 
this range [24| • A normalization option has to be chosen 
to ensure agreement with the pp experimental data and 
this typically leads to reductions of about 15% at our two 
lower energies but almost 20% at 2.83 GeV. To eliminate 
any ambiguities, the values that we have assumed for the 
pp elastic cross section, integrated over the solid angle 
defined by 5.0° < 6*i a b < 8.5°, are presented in Table U 
for our three energies. 

TABLE I: Total cross sections for the pp — > pp{K + K^}^^.^ 
and pp— >ppij> reactions at three beam energies T p and corre- 
sponding excess energies ekk with respect to the ppK + K~ 
threshold. The excess energy with respect to the <f> threshold 
is given by e$ = ekk — 32.1 MeV. The cf> cross section has 
been corrected for the 49.1% branching ratio 4]. For both 
channels the first error is statistical and the second system- 
atic. The values that we have taken for a pp , the integral of the 
proton-proton elastic cross section over the solid angle range 
5.0° < 6*i a b < 8.5°, are also given. The associated w ±6% un- 
certainty has not, however, been compounded with the other 
errors. 



T P 


e K+K- 




O" llon -0(tOt) 


(70 (tot) 


[GeV] 


[MeV] 


[mb] 


[nb] 


[nb] 


2.65 


51 


5.52 


16±1±1 


33 ± 2 ± 4 


2.70 


67 


5.44 


30 ± 2 ± 3 


64 ± 4 ± 10 


2.83 


108 


5.06 


98 ± 8 ± 15 


133 ±12 ±27 



The value that we obtain for the luminosity at 
2.65 GeV on the basis of the pp cross section given in 
the table could be checked using a novel technique [25| . 
This depends upon the energy loss of the coasting pro- 
ton beam caused by its multiple traversals through the 
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FIG. 10: Differential cross sections for the pp — > ppK + K~ reaction at 2.65 GeV with respect to the invariant masses of the 
K + pp and the K~pp systems and their ratio. The data are divided (a) into a 0-poor region with M(K + K~) < 1.01 GeV/c 2 
(upper panels) and (b) rich with M(K + K~) > 1.01 GeV/c 2 (lower panels). The predictions obtained using the standard 
parameters are shown with no fsi effect for the K~ from the decay of the <f> (solid curve) and 20% (dashed). The simulations 
that describe the two-particle ratios also represent well these three-body ratios. 



very thin target. The resultant frequency change can 
be measured with high accuracy using the spectrum of 
the Schottky noise. Combined with measurements of the 
beam current, this yielded luminosities with an expected 
precision of about ±6%. The two methods give consis- 
tent results within this uncertainty and, although the 
frequency-change technique could not be applied as reli- 
ably at the higher energies, it is expected that the overall 
luminosity uncertainty is on about this level. 



Our total cross section results for non-0 production 
are plotted in Fig. [T2] along with results taken from 
DISTO ll) and COSY-11 0, HE S3- It is seen that 
four-body phase space cannot describe simultaneously 
the pp — ► ppK + K~ data at high and low excess energies. 
The situation is improved only a little if account is taken 
of the strong attraction between the two protons in the 
^So state but a greater improvement is achieved through 
the introduction of the fsi in the K~p system. Never- 
theless, the COSY-11 data at excess energies of 20MeV 
or less still seem to be underestimated. This might be 
connected to the incomplete description of the K + K~ 
spectra of Fig. [3] at very low invariant masses. However, 
it must be stressed that in the evaluation of the COSY-11 
acceptance no account was taken of any K~p fsi. 



VIII. CONCLUSIONS AND OUTLOOK 

We have investigated the pp — ► ppK + K~ reaction at 
three energies. In those parts of the K + K~ mass spectra 
where <j> production is small, we find a marked difference 
between the K~p and K + p invariant-mass distributions. 
There is a strong peaking in the K~~p/K + p ratio towards 
the Kp threshold for both the COSY-11 and our data. 
This clearly indicates that the pp — > ppK + K~ reaction 
cannot be dominated by the undistorted production of 
a single scalar resonance ao or /o- Either a mixture of 
the two is produced or one of the kaons from the ao / /o 
decay interacts with a final proton. As we have shown, 
these data can be explained quantitatively through the 
introduction of a simple K~p final state interaction. 

In kinematic regions where <fi production is significant, 
the K~p/ K + p asymmetry is reduced because the meson 
may travel some distance before decaying. Our data in 
such conditions are well fit by assuming that about 20% 
of the K~ from the <j) decay interact with one or both 
final protons. The uncertainty in this is, however, quite 
large because of the limited precision with which the cf> 
contribution to the K + K~ mass spectrum can be iso- 
lated. 

There is no universal technique for modelling the fsi 
when more than one pair of particles is involved. The 



10 



-a 

73 



10" 1 



(a) 



■4t 




(b) 




2.4 



2.42 




u 2 4 



:.t 2.42 2.44 2.46 
K p p invariant mass [GeV/c 2 ] 

FIG. 11: The ratio of the differential cross sections of the 
pp — > ppK + K~ reaction with respect to the K~pp to K + pp 
invariant masses for (a) the <^-poor region and (b) the 0-rich 
at 2.70 GeV with (c) and (d) showing respectively the analo- 
gous results at 2.83 GeV. The predictions obtained using the 
standard parameters are shown with no fsi effect for the K~ 
from the decay of the <j> (solid curve) and 20% (dashed) . The 
simulations that describe the two-particle ratios also repre- 
sent well these three-body ratios. 



product ansatz has the benefit that the same simula- 
tion with the same parameters reproduces also the strong 
peaking in the ratio of the K~pp/K + pp mass distribu- 
tions. 

We find an effective K~p scattering length that is of 
the same order of magnitude as that needed to describe 

F 1 

the free scattering data [19(. However, the error bars are 




[MeV] 



FIG. 12: Energy dependence of the non-<^> contribution to 
the pp — > ppK + K~ total cross section. In addition to our 
three points (closed circles), there is a high energy point from 
DISTO (open circle) [lj and four points close to threshold 
from COSY-11 (closed squares) (open triangle) [2(|, and 
(open square) [23] • The dashed curve represents the energy 
dependence from four-body phase space. The dot-dashed in- 
cludes the effects of the pp final state interaction whereas the 
solid curve contains in addition distortion from the K~p fsi, 
as described in the text. 



large and our data are primarily sensitive to \a\ rather 
than its phase. It is important to stress that we do not 
know the relative strengths for the production of isospin 
1 = and 1 = 1 K~p pairs in the pp — > ppK + K~ reac- 
tion and hence the relative weights of ao and a\. Some 
information on this might come from looking at the data 
on the pp -> AT+pE°(1385) and pp -> #+pA(1405) re- 
actions taken at 2.83 GeV [g]. However, even this would 
not give access to the phase between the 1 = and 1=1 
amplitudes. It must also be stressed that we have ne- 
glected any energy dependence of the parameter a. An- 
other important caveat is that we have not considered 
any possible interaction of the K + with the protons. Al- 
though there is no evidence for any strong effect of this 
nature in the pp — > K + Ap reaction [28| . since the K + p 
interaction is repulsive, its neglect might be interpreted 
as extra attraction in the K~p system. 

We have clearly shown the importance of the 
K~p interaction, which is possibly related to the A(1405) 
or the E°(1385), but we must now ask whether there is 
any evidence for effects that might be connected with 
the production of the ao//o resonances. The simula- 
tion of the energy dependence of the total cross sec- 
tion for the production of the non-<^> component in the 
pp — > ppK + K~ reaction underestimates the COSY-11 
data very near threshold. Even more intriguing, and 
possibly connected, are the lowest mass points in the 
K + K~ spectra. At all three energies shown in Fig. [3l 
these points lie much higher than the simulations and a 
similar behavior is visible in the DISTO results [l|as well 
as in our data on the pn — > dK + K~ reaction 31. The 
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mass scale of the variation is not that of the widths of the 
scalar resonances, which are quite large. It is tempting to 
suggest that this structure might be due to the opening 
of the K°K° channel at a mass of w 8MeV/c 2 , which 
induces some cusp structure that changes the energy de- 
pendence of the total cross section near threshold 0, H^] • 
This would require a very strong K + K~ =± K°K° chan- 
nel coupling, which might be driven by the ao//o res- 
onances. Thus, although the pp — > ppK + K~ reaction 
may not be ideal for investigating the properties of scalar 
states, their indirect effects might still be crucial. 

Given that the K~ is strongly attracted to protons, 
it is natural to expect that similar fsi effects should 
exist for other reactions. The K°d/K + d ratio in the 
pp — » K + K°d reaction does indeed show a very strong 
enhancement near the Kd threshold Q- The interpre- 
tation is, however, slightly more complex because spin- 
parity constraints means that there must be at least one 
p-wave in the final state, even at low energies. Unfortu- 
nately, the MOMO data onpd -> K + K~ 3 He [H cannot 
be used to isolate any K~ 3 He interaction because the 
charges of the kaons were not identified and only aver- 
aged spectra could be studied [3ll ]. 

There is currently an overall uncertainty in our val- 
ues of the cross sections because of the limitation in the 
proton-proton elastic scattering database in our region 
of energy and angle, though we hope that this will be al- 
leviated by future measurements [25| . With our new esti- 
mates for the luminosity, the ^-production cross sections 
are reduced compared to those of our first analysis Q, 
which relied upon the non-renormalized 2004 SAID so- 
lution [23j . Furthermore, the better separation of the (j> 
and non-<^> events goes by chance in the same direction. 
These two effects lower the ratio of (f> to uj production in 
proton-proton collisions near threshold so that it is now 
only about a factor of six above the Okubo Zwcig Iizuka 
limit [13]. 
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APPENDIX A: FINAL STATE INTERACTIONS 

When two particles interact strongly in the final state, 
the resulting matrix element involves an average of a pro- 
duction operator with the relative wave function ip$(r) of 
the strongly interacting pair. If the interaction is of very 
short range, the wave function may be evaluated at the 
origin to leave an enhancement factor 

F 2 (g)«V*(0)cx^, (Al) 

where <f is the relative momentum in the pair and D(q) 
is the 5-wave Jost function 18]. In the commonly used 
scattering length approximation, one retains only the lin- 
ear term in q, in which case F(q) = 1/(1 — iqa), where a 
is the scattering length of the interacting pair. 

On the other hand, if two or more pairs of particles 
interact in the final state, there is no reliable prescrip- 
tion to evaluate an analogous enhancement factor since 
a three-body equation then needs to be solved. If we de- 
note the interacting pairs as 12 and 13, the corresponding 
wave function will be ^$ 12 rf 13 {ri2, ^13)- We now make the 
ad hoc assumption that this wave function factorizes in 
the form 

*?12 ,913 (712 ,f 13 ) = 1pq 12 (fi 2 ) X ijjfa (fi 3 ) . ( A2) 

In this case, the three-body enhancement factor is simply 
the product of the two-body factors, evaluated at the 
appropriate relative momenta: 

-^3(912, 913) = F 2 {qi 2 ) x F 2 {qi 3 ) 



(1 - i<7i 2 ai2)(l - %3ai3) ' 

It must be stressed that Eq. (|A3|) is merely an ansatz 
to try to understand our ensemble of data and, even if 
it provides a satisfactory description of these, this does 
not mean that it is applicable more generally. However, 
it has also been used to treat the pp — > pprj reaction near 
threshold, where all the final pairs of particles interact 
strongly (l7| . 

The approach does retain the necessary pole structure 
when q\2 = —i/a\2 and similarly for q\ 3 . Thus particle- 
1 can interact simultaneously with both 2 and 3. Fur- 
thermore, to lowest order in the momenta, the ansatz 
corresponds to the scattering length of particle-1 from 
a composite 2+3 with the desired combined scattering 
length of a = a\ 2 + ai3- 

Although the S'-wave pp interaction is well known, 
the K~p interaction is far more complex because of the 
channel couplings to S7r and also because there are two 
isospins 1 = and 1=1. In a recent study, within 
a chiral SU(3) unitary approach [Hit]. K~p data leading 
to a variety of channels have been fit and values of the 
(complex) K ~p scattering lengths deduced for 1 = and 
1 = 1. Their full model resulted in 

a = (-1.64 + i0.75)fm, 

ai = (-0.06 + z0.57)fm. (A4) 
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These values describe the rapid drop in the K~ p elastic 
amplitudes from threshold. 

If isospin breaking induced by the significant mass dif- 
ferences is neglected, the K~p scattering length a^-p = 
(<jq + a±)/2. However, this is not necessarily the pa- 
rameter that is relevant for the analysis of the pp — > 



pK + {K~p} because it is not clear whether the reac- 
tion mechanism preferentially excites I = or I = 1 
{K~p} states or some combination thereof. The values in 
Eq. (|A4[) should therefore be considered merely as order- 
of-magnitude estimates of what might be expected when 
making fits to our data. 
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